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Proton Radiation Effects and Metrics 


Total Ionizing Dose (TID) - cumulative damage 
resulting from ionization (electro n-hole pair^ 
formation) causing 

- Threshold voltage shifts 

Metric used for TID studies: 

Dose = energy deposited per uniHnas^^^^ 
material that comprises sensitive volume 

- 1 Gy = 1 J/kg 

- 1 rad = 100 erg/g 
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Proton Radiation Effects and Metrics 


Displacement Damage - cumulative damage resulting from 
displacement of atoms in semiconductor lattice structure 
causing: SSS 

- Carrier lifetime shortening °o 0 o^o 0 o 0 c 

- Mobility degradation 0 o^S°o 0 o 0 c 

- Charge transfer degradation ^ . % o c°c°c 0 ( 

in imaging devices 

Two metrics used for displacement damage studies: 

- Displacement Damage Dose = energy going into displaced 
atoms (nonionizing energy) per unit mass of material that 
comprises sensitive volume 

- Equivalent Proton Fluences: 10 MeV often-used standard; 
other energies can be used 


• interstitial 
O Vacancy 
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Proton Radiation Effects and Metrics(^ 

* Single Event Effects (SEE) - event caused by 
single incident proton 

- Non-destructive - SEU, SET, MBU, SHE 

- Destructive - SEL, SEGR, SEB 

* Results commonly presented as function of 
incident proton energy 

- Most SEE result from nuclear reaction products, which 


are influenced significantly by incident proton energy. 



GCR Radiation Effects and Metrics 



• Single Event Effects (SEE) - event caused by single 
incident ion 


- May be caused by direct ionization (usually the case for 
incident heavy ions) 


- May be caused by nuclear reaction products (usually the case 
for incident protons) 



Metric commonly used for heavy ion induced SEE is Linear 
Energy Transfer (LET) 

LET = energy lost by ionizing particle per unit path length in 
sensitive volume 

- Path length often expressed as areal density by dividing by 

material density ■ 

- LET units commonly used are MeV-cm 2 /mg. 

- For some situations, especially small feature .^fe.ffatg^ ^ 
size devices, LET may not be sufficient. 
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Electron Radiation Effects and Metrics 


TID - similar to that for protons 

Displacement Damage 

- Do less damage than equal number of protons 

- Metrics similar; 1 MeV equivalent electron fluences are 
often used. 

Charging/Discharging Effects 

- Surface charging caused mainly by low B 

- energy electrons/plasma 

- Deep dielectric charging caused mainly by high energy 
electrons 

- Key parameter is potential difference induced by 
charging between dielectric and conductive surface. 


Introduction 


The radiation environment must be 
understood and accurately modeled 

- reliable, cost-effective designs 

- implement new space technologies 
Underestimating radiation levels leads to 

- excessive risk 

- degraded system performance 

- loss of mission lifetime 
Overestimating radiation levels leads to 

- excessive shielding 

- reduced payloads 

- over-design 

- increased cost 
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The Solar Activity Cycle 




The sun is a source and modulator 
of space radiations. 

Its approximately 11 -year cycle 
typically consists of 

- 7 years of solar maximum when 
activity levels are high 

- 4 years of solar minimum when 
activity levels are low 

Sunspot numbers are commonly 
used indicators of solar activity. 




Galactic Cosmic Rays 


• Galactic Cosmic Rays (GCR) are high-energy 
charged particles that originate outside our solar 
system that are present at low flux levels. 
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GCR Properties 


Composed mainly of 
hadrons 

Consist of all naturally 
occurring elements* 

- 87% protons 

- 12% alpha particles 

- 1% heavier ions 
Energies: up to 10 11 GeV! 
Fluxes: 1 to 10 cm^s 1 


The relative abundance of 
i Galactic Cosmic Rays - 


Variation with Solar Cycle 


-Solar Maximim 
-Solar Minimum 


X Protons 

^ \ 

Helium \ 


- Solar Minimum 

- Solar Maximum 
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GCR Models 


NASA and MSU models originated independently 

- Both based on theory of solar modulation 

- Describes penetration of GCR into heliosphere from outside 
and transport to near earth 

- Solar modulation results in variation of GCR fluxes over solar 
cycle 

Implementation of solar modulation differs 

- NASA model determines solar modulation from near earth 
GCR measurements, including detection of secondary 
neutrons with ground-based monitors. 
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Solar Particle Events 



• Solar flares 


- Occur when localized energy storage in coronal 

magnetic field becomes too great and burst of energy is 
released 



Coronal Mass Ejections (CMEs) 

- Large eruption of plasma that drives a shock wave 
outward and accelerates particles 

Responsible for major disturbances 

- Interplanetary space 

- In earth’s magnetosphere I * 


or* CempoiswiisSndi Syst 


Characteristics of CMEs 

• Hadron composition* 
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- 96.4% protons 

- 3.5% alpha particles 

- 0.1% heavier ions (not to be neglected!) 

• Energies: up to ~ GeV/nucleon 

• Event magnitudes: 

- > 10 MeV/nudeon integral fluence: can exceed 10 9 cm* 2 

- > 10 MeV/nudeon peak flux: can exceed 10 s cmV 

* D.V. Reames, Space Sci. Rev., 1999 


To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components and Systems 
(RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 


9 







Probabilistic phenomena 
Truncated power law 
function describes 
essential features of 
distribution of event 
fluences: 

- Sudden bursts of energy 
released that span orders 
of magnitude 

- Smaller event sizes follow 
power law function 


D 7 10 8 1 

> 30 MeV Fluence (crrv 2 ) 


ii mil!. 
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Cumulative Fluence Models 

Solar Protons 


• Once the distribution of 
event fluences is known 
the cumulative fluence 
during a mission can be 
calculated. 

• JPL91 model 

- Uses Monte Carlo 
simulations 

• Other Monte Carlo 
variations 

- ONERA, Nymmik 

• Analytic 
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Cumulative Fluence Models 

Solar Heavy Ions 


Preliminary model by Tylka 
for 2 energy bins each of He, 
CNO group and Fe 
PSYCHIC model of NASA 
GSFC 

- Statistical model of alpha 
particles based on 28 years 
of data from IMP-8 and GOES 

- Major heavy ions C through 
Fe determined from ACE SIS 
instrument measurements. 

- Remaining minor elements 
scaled according to ISEE-3 
measurements and corrected 
photospheric abundance 
model 
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Worst Case Events 

CREME96 Model based on October 1989 event 
and provides 3 intensity levels 
- Peak flux, worst day, worst week 
Comparisons to QinetiQ’s CREDO experiment 
onboard MPTB show that 3 events during solar 
cycle 23 approximately equaled the “worst day” 
model _ 8 


. : ^ 

| 1Xl ° 2 

3 

't 1x10' 
I? 1x10° 

, C 

. . . ' > . IxlO 1 

? If: m % h $ ' 1 



1 2 ~+- 2926 (start) 
. (2 3 -*r 2927 

* ‘X. 

8 2929 (peak) 

5 - ♦ 2931 

*v\ 

3 , A 2932 

i 2933 


CREME96 Worst Day 






fills hi in 


' ■ V ■ s 

: I -.-f' f •.-!■ , i ^ | 


i ■ i " ' I ■ I 




1 


Solar and Galactic Cosmic Ray 

Space Environment ** 


• Three situations that need 
to be considered for 
spacecraft design and for 
SEE rate predictions: 

- Ambient Galactic Cosmic 
Ray flux 

- Cumulative solar particle 
event fluence 

- Worst case solar particle 
event flux 
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Future Challenges 
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• Strategies for manned and robotic missions for 
new interplanetary exploration initiatives 

- How to address our limited knowledge of future solar 
cycle activity. Solar activity has a direct influence on 
GCR levels, which are difficult to shield against 

■ Astronaut exposure 

■ SEE in high performance microelectronics 


Future Challenges 



• Addressing the lack of 

predictability of solar particle 
events for new manned 
exploration initiatives requires 
more sophisticated strategy 
- Establish measurement system 
in inner heliosphere for early 
detection of events. 



- Subsequent prediction of 
properties such as arrival time, 
duration, intensity and energy 
spectrum must be made for 
missions to the moon, Mars and 
possibly beyond. 
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1 Introduction 

Before the space era the only manifestations of the presence of radiations in space were 
the deformation of the ionized tail of comets caused by the solar wind, the aurora borealis whose 
origin was not well understood and the ionization of air, secondary cosmic-ray showers and the 
isotopes (carbon 14 for example) produced by cosmic radiation. We had to wait until the 
beginning of the space conquest to discover in 1958 the presence around the Earth of very-high 
energy charged particles (Van Allen belts). Since then, it has become evident that the space 
environment is a highly aggressive medium. Beyond the natural protection provided by the 
Earth's atmosphere, various types of radiation can be encountered. Their characteristics (energy 
and nature), origins and distributions in space are extremely variable. This environment degrades 
electronic systems and on-board equipment in particular and creates radiobiological hazards 
during manned space flights. 

Based on several tens of years of the space adventure, a detailed analysis of the problems 
on satellites shows that the part due to the space environment is not negligible. It appears that the 
malfunctions are due to problems linked to the space environment (9 to 21%), electronic 
problems (6 to 16%), design problems (11 to 25%), quality problems (1 to 8%), other problems 
(1 1 to 33%) and problems that are still unexplained (19 to 53%) [Space, 1994], It is clear that the 
unexplained problems are either problems linked to the space environment, or to the electronics, 
or to the design, or to the quality or otherwise but the information collected on the ground is 
generally/often not sufficient to define the origin of the problem. The space environment is 
largely responsible for about 20% of the anomalies occurring on satellites and a better 
knowledge of that environment could only increase the average lifetime of space vehicles. 

This naturally leads to a detailed study of the space environment and of the effects that it 
induces on space vehicles and astronauts. The nature of this environment varies greatly between 
low orbits and the higher altitudes such as the geostationary orbit and beyond. Among its 
components, we only examine the ionizing charged particles here, i.e. the particles trapped in the 
radiation belts, the solar flare and cosmic radiation particles. From the point of view of the 
effects, the degradations will differ according to the energy of the particles, to their nature and to 
the satellite orbit. The degradations and disturbances induced by space radiation in the materials 
and the electronic components are phenomena that have been studied for many years. Two 
categories of effects should be noted: 

- the cumulative effects such as the aging of thermal control coatings, of the optics and 
electronics and the erosion of materials; 

- the sporadic effects such as noises in the detectors and optics, single event effects in 
highly integrated electronic circuits and electrostatic discharges. 

The increasingly frequent presence of man in space and the projects such as the distant 
and long-duration missions (Lunar base, flight to Mars, etc.) pose the problem of the biological 
effects induced, essentially in the long term, by high energy radiation. Radiation can have two 
possible types of biological effects: 

- immediate, permanent or delayed non stochastic effects (destruction or modification of 
cells), the speed with which the symptoms appear and their seriousness increase in proportion to 
the exposure to the radiation; 

- stochastic effects, associated with the modifications to the cells, whose probability of 
appearing in the long term increases in proportion to the irradiation (cancers, leukemia, genetic 
effects). 
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2 The Solar Activity Cycle 

The sun is both a source and a modulator of space radiations. Understanding its cyclical 
activity is an important aspect of modeling the space radiation environment. The solar activity 
cycle is approximately 1 1 years long. During this period there are typically 7 years during solar 
maximum when activity levels are high and 4 years during solar minimum when activity levels 
are low. In reality the transition between solar maximum and solar minimum is a continuous one 
but it is often considered to be abrupt for convenience. At the end of each 11 -year cycle the 
magnetic polarity of the sun reverses and another 1 1 -year cycle follows. Thus, strictly speaking 
the total activity cycle is approximately 22 years long. Of the space radiations considered here 
the magnetic polarity apparently only affects the galactic cosmic ray fluxes [Ba96a], and not the 
trapped particle or solar particle event fluxes. Thus, things are often viewed on an approximately 
1 1 -year cyclical basis. 

Two common indicators of this approximately 1 1 -year periodic solar activity are sunspot 
numbers and solar 10.7 cm radio flux (F 10 . 7 ). The most extensive record is that of observed 
sunspot numbers, which dates back to the 1600s. This record is shown in Figure 1. The 
numbering of sunspot cycles began in 1749 and it is currently near the end of solar cycle 23. The 
record of F\q j began part way through solar cycle 18 in the year 1947 and is shown in Figure 2. 


Yearly Averaged Sunspot Numbers 1610-1998 
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Year 

Figure 2. Measured values of solar 10.7 cm radio flux. 

Although sunspot numbers and F 10.7 are commonly accepted indicators of solar activity, 
quantitative relations to measured radiation events and fluxes are not necessarily straight 
forward. Large solar particle events are known to occur with greater frequency during the 
declining phase of solar maximum [Sh95]. Trapped electron fluxes also tend to be higher during 
the declining phase [Bo03]. Trapped proton fluxes in low earth orbit (LEO) reach their 
maximum during solar minimum but exactly when this peak is reached depends on the particular 
location [Hu98]. Galactic cosmic ray fluxes are also at a maximum during solar minimum but in 
addition depend on the magnetic polarity of the sun [Ba96a]. 

There has been considerable effort put into forecasting long-term solar cycle activity. A 
review of a number of the methods is presented by Hathaway [Ha99], These include regression 
methods, which involve fitting a function to the data as the cycle develops. Also discussed are 
precursor methods, which estimate the amplitude of the next cycle based on some type of 
correlation with prior information. These methods can also be combined. In addition, physically 
based methods are being developed based on the structure of the magnetic field within the sun 
and heliosphere [Sc96], [Di06], 

However, accurate methods for predicting future solar cycle activity levels prior to the 
start of the cycle have thus far been elusive. A potential breakthrough, however, has recently 
been reported that uses a combination of computer simulation and observations of the solar 
interior from instrumentation onboard the Solar and Heliospheric Observatory (SOHO) [Di06]. 
Given the current state of this modeling, probabilistic models of solar activity can be useful. 
Such a model of F 10.7 is shown in Figure 3 [Xa02]. This also illustrates the general behavior of 
the observed cyclical properties, at least over recent cycles. The greater the peak activity of a 
cycle, the faster the rise-time to the peak level. Furthermore the cyclical activity is asymmetric 
such that the descending phase of the cycle is longer than the ascending phase. 


1-5 




To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components 
and Systems (RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 


6 



Time from Start of Solar Cycle (Years) 

Figure 3. Probabilistic model of F 10 . 7 . The various curves are labeled as a function of 
confidence level that the activity shown will not be exceeded [Xa02]. 


3 Galactic Cosmic Rays 
3.1 General Characteristics 

Galactic cosmic rays (GCR) are high-energy charged particles that originate outside of 
our solar system and are believed to be mainly remnants from supernova explosions. These 
explosions occur about once every 50 years in our galaxy. Some general characteristics of GCR 
are listed in Table 1 . They are composed mainly of hadrons, the abundances of which are listed 
in the Table. A more detailed look at the relative abundances is shown in Figure 4. All naturally 
occurring elements in the Periodic Table (up through uranium) are present in GCR, although 
there is a steep drop-off for atomic numbers higher than iron (Z=26). Energies can be as high as 
10 11 GeV, although the acceleration mechanisms to reach such high energies are not understood. 
Fluxes are generally a few cm'V 1 , and vary with the solar cycle. Typical GCR energy spectra 
for a few of the major elements during solar maximum and solar minimum are shown in Figure 
5. It is seen the spectra tend to peak around 1 GeV per nucleon. The flux of the ions with 
energies less than about 10 GeV per nucleon is modulated by the magnetic field in the sun and 
solar wind. During the high activity solar maximum period there is significantly more 
attenuation of the flux, resulting in the spectral shapes shown in Figure 5. 
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Hadron 

Composition 

Energies 

Flux 

Radiation Effects 

Metric 

87% protons 
12% alphas 
1% heavier ions 

Up to -10“ GeV 

1 to 10 cm' 2 s' 1 

SEE 

LET 


Table 1 . Characteristics of Galactic Cosmic Rays. 



Atomic Number (Z) 

Figure 4. Abundances of GCR up through Z = 28. 
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Energy (MeV/amu) 

Figure 5. GCR energy spectra for protons, helium, oxygen and iron during solar maximum and 

solar minimum conditions [Ba96a]. 


Single Event Effects (SEE) are the main radiation effects caused by GCR in 
microelectronics and photonics. The metric traditionally used to describe heavy ion induced 
SEE is linear energy transfer (LET). LET is the energy lost by the ionizing particle per unit path 
length in the sensitive volume. For SEE studies the path length is often divided by the material 
density and expressed as an areal density. The units of LET that are commonly used are then 
MeV-cm 2 /mg. 

For SEE analyses energy spectra such as those shown in Figure 5 can be converted to 
LET spectra. Such integral LET spectra for solar maximum and solar minimum conditions are 
shown in Figure 6. These spectra include all elements from protons up through uranium. The 
ordinate gives the flux of particles that have an LET greater than the corresponding value shown 
on the abscissa. Given the dimensions of the sensitive volume this allows the flux of particles 
that deposit a given amount of charge or greater to be calculated in a simple approximation. In 
some modem devices, the LET metric may not be useful because of the highly scaled nature of 
devices, the complexity of the geometry, or the increased use of metal over-layers. In addition, if 
nuclear reactions play a significant role in producing SEE, the LET metric is not valid for this 
situation. Exposure of astronauts to GCR is a serious consideration for manned missions 
because GCR are difficult to shield against given the typical weight constraints of missions. 


1-8 



To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components 
and Systems (RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 


9 



The LET spectra shown in Figure 6 are applicable to geosynchronous and interplanetary 
missions where there is no geomagnetic attenuation. The earth’s magnetic field, however, 
provides significant protection. Due to the basic interaction of charged particles with a magnetic 
field, the charged particles tend to follow the geomagnetic field lines. Near the equator the field 
lines tend to be parallel to the earth’s surface. Thus all but the most energetic ions are deflected 
away. In the polar regions the field lines tend to point toward the earth’s surface, which allows 
much deeper penetration of the incident ions. The effect of the geomagnetic field on the incident 
GCR LET spectrum during solar minimum is discussed for various orbits in [Ba97]. 

3.2 Galactic Cosmic Ray Models 

The original Cosmic Ray Effects in MicroElectronics (CREME) suite of programs of 
Adams [Ad87] was developed specifically for microelectronics applications. It turned out to be 
a very useful and popular tool and has been updated since then. CREME96 is the current version 
[Ty97] and uses the GCR model of Moscow State University (MSU) [Ny96a], Recent work has 
begun on a subsequent update [Ad07]. 

In principle the MSU model is similar in approach to a GCR model that was originated 
independently at NASA by Badhwar and O’Neill [Ba96a]. Both models are based on the 
diffusion-convection theory of solar modulation [Pa85]. This is used to describe the penetration 
of cosmic rays into the heliosphere from outside and their transport to near earth at 1 
Astronomical Unit (AU). The solar modulation is used as a basis to describe the variation of 
GCR energy spectra over the solar cycle, as shown in Figure 5. However, the implementation of 
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the solar modulation theory for the two models is different. The Badhwar and O’Neill model 
estimates the modulation level from GCR measurements at 1 AU. Correlations to ground-based 
neutron monitor counting rates are then made to establish long-term predictive capability. The 
MSU model is not as direct but uses multi-parameter fits to ultimately relate solar cycle 
variations in GCR intensity to observed sunspot numbers. 

Comparisons of the GCR proton and alpha particle spectra of the two models above plus 
that used in the QinetiQ Atmospheric Radiation Model (QARM) show discrepancies among all 
three models for narrow time ranges [Le06]. Examples of this are shown in Figure 7 for protons. 
This is not surprising considering the details of the solar modulation implementation are 
different. However, similar predictions are seen for the total fluence over the course of a solar 
cycle. 




Figure 7. GCR proton energy spectra predicted by the MSU, Badhwar and O’Neill, and QARM 

models for two different dates [Le06]. 


The recent high-quality measurements of GCR heavy ion energy spectra taken on the 
Advanced Composition Explorer (ACE) satellite make possible an interesting test of the GCR 
models. Comparisons of model results and the ACE data for the 1997 solar minimum period are 
shown in Figure 8 for 4 of the major elements in the energy range of about 50 to a few hundred 
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MeV per nucleon. The NASA results incorporate a recent update [ON07]. It is seen that both 
models yield good results for heavy ions although the updated NASA model is more accurate for 
this situation in terms of spectral shape and root-mean-square deviation from the data. 



Figure 8. Comparison of the updated NASA model of Badhwar and O’Neill and the MSU model 
to measurements made with instrumentation onboard the ACE satellite during 1997. After 

[DaOl], 


A recent development led by the California Institute of Technology is to use a transport 
model of GCR through the galaxy preceding the penetration and subsequent transport in the 
heliosphere. [DaOl]. During the initial propagation of GCR through the galaxy use is made of 
knowledge of astrophysical processes that determine the composition and energy spectra of 
GCR. Comparisons of the fitted model spectra to the ACE satellite measurements are shown in 
Figure 9. The model spectra do not reflect solar modulation so the model is still a work in 
progress. The fitting is done to demonstrate that the model has the potential to closely reproduce 
the ACE measurements. The elements C and Fe are GCR primaries while B, Sc, Ti and V are 
GCR secondaries produced by fragmentation of primaries on interstellar H and He. The goal of 
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this new approach is to provide an improved description of GCR composition and energy spectra 
throughout the solar cycle. 



Figure 9. Demonstration that the new approach of the California Institute of Technology can 
describe GCR energy spectra measured by ACE instrumentation during 1997. The model does 

not yet incorporate any solar modulation [DaOl]. 


4 Solar Particle Events 

4.1 General Characteristics 

It is believed that there are 2 categories of solar particle events and that each one 
accelerates particles in a distinct manner. Solar flares result when the localized energy storage in 
the coronal magnetic field becomes too great and causes a burst of energy to be released. They 
tend to be electron rich, last for hours, and have an unusually high 3 He content relative to 4 He. A 
Coronal Mass Ejection (CME), on the other hand, is a large eruption of plasma (a gas of free ions 
and electrons) that drives a shock wave outward and accelerates particles. CMEs tend to be 
proton rich, last for days, and have a small 3 He content relative to 4 He. A review article by 
Reames gives a detailed account of the many observed differences between solar flares and 
CMEs [Re99], 
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CMEs are the type of solar particle events that are responsible for the major disturbances 
in interplanetary space and the major geomagnetic disturbances at earth when they impact the 
magnetosphere. The total mass of ejected plasma in a CME is generally around 10 15 to 10 17 
grams. Its speeds can vary from about 50 to 1200 km/s with an average speed of around 400 
km/s. It can take anywhere from about 12 hours to a few days to reach the earth. Table 2 lists 
some further general characteristics of CMEs. 


Hadron 

Composition 

Energies 

Integral Fluence 
(>1 OMeV/nucleon) 

Peak Flux 
(>1 OMeV/nucleon) 

Radiation 

Effects 

96.4% protons 




TID 

3.5% alphas 

Up to 

Up to ~10 10 cm' 2 

Up to ~10 5 cm'V 1 

DD 

~0.1% heavier 

~GeV/nucleon 



SEE 

ions 






Table 2. Characteristics of CMEs. 


All naturally occurring chemical elements ranging from protons to uranium are present in 
solar particle events. They can cause permanent damage such as Total Ionizing Dose (TID) and 
Displacement Damage (DD) that is due mainly to the proton and possibly the alpha particle 
component. Just because the heavy ion content is a small percentage does not mean it can be 
ignored. Heavy ions, as well as protons and alpha particles in solar particle events, can cause 
both transient and permanent SEE. 

Figures 10 and 11 illustrate the periodic yet statistical nature of solar particle events. 
They are plots of the daily solar proton fluences measured by the Interplanetary Monitoring 
Platform-8 (IMP-8) and Geostationary Operational Environment Satellites (GOES) over an 
approximately 28 year period. Figure 10 shows > 0.88 MeV fluences while Figure 11 shows > 
92.5 MeV fluences. The solar maximum and solar minimum time periods are shown in the 
figures to illustrate the dependence on solar cycle. 
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Figure 10. Daily fluences of > 0.88 MeV protons due to solar particle events between 

approximately 1974 and 2002. 



Figure 11. Daily fluences of > 92.5 MeV protons due to solar particle events between 

approximately 1974 and 2002. 
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The available solar particle data that cover the largest period of time are for protons. Since 
the available solar heavy ion data are not nearly as extensive, solar proton models and solar 
heavy ion models will be discussed separately. 

4.2 Solar Proton Models 

Sections 4.2. 1 - 4.2.5 describe the application of probabilistic methods to solar proton 
event data, including their origin. Section 4.2.1 establishes the probabilistic nature of events. 
Section 4.2.2 then describes the distribution of event magnitudes. 4.2.3 and 4.2.4 describe 
modeling cumulative fluences over the course of a mission. 4.2.5 discusses worst-case events 
during a mission. 

4.2.1 Self-Organized Criticality and the Probabilistic Nature of the Energy Release 

Process 

Substantial efforts have been put into studies of the occurrence of solar particle events. 
One of the main goals is to find a reliable predictor of events. Despite this significant 
international effort, solar particle events can occur suddenly and without obvious warning. In 
addition to potential problems with electronic systems and instrumentation, this is an especially 
serious concern for new space initiatives that plan to send manned spacecraft to the moon, Mars 
or interplanetary space. Thus, there is strong motivation to develop predictive methods for solar 
particle events. It is hoped that the apparent stochastic character can be overcome and 
predictability achieved if precursor phenomena such as x-ray flares or magnetic topology 
signatures can be properly interpreted or if the underlying mechanisms are identified. Whether 
the nature of the energy release process for solar particle events is deterministic or stochastic is a 
very basic question. More specifically, the question is whether it is possible to predict the time 
of occurrence and magnitude of solar particle events or if probabilistic methods are necessary. 

The self-organized criticality (SOC) model is a phenomenological model originated by 
Bak, Tang and Wisenfeld [Ba87] that can give insight into the basic nature of a system. It 
postulates that a slow continuous build-up of energy in a large interactive system causes the 
system to evolve to a critical state. A minor, localized disturbance can then start an energy- 
releasing chain reaction. Chain reactions and therefore energy releasing events of all sizes are an 
integral part of the dynamics, leading to a “scale invariant” property for event sizes. This scale 
invariance results in power function distributions for the density functions of event magnitudes 
and waiting times between events. As a result of this basic nature it is generally assumed in the 
literature that accurate predictions of the magnitude and time of occurrence of such events are 
not possible. A system in a SOC state is therefore generally assumed to be probabilistic in 
nature. 

Applications for the theory of SOC have been found in natural phenomena such as 
earthquakes, avalanches and rainfall. A useful conceptual aid is the sandpile. If sand is dropped 
one grain at a time to form a pile, the pile soon becomes large enough that grains may slide down 
it, thus releasing energy. Eventually the slope of the pile is steep enough that the amount of sand 
added is balanced, on average, by the amount that slides down the pile. The system is then in the 
critical state. As single grains of sand are subsequently added, a broad range of consequences is 
possible. Nothing may happen or an avalanche of any size up to a “catastrophic” one may occur. 
The dynamics of this interactive system do not allow accurate predictions of when an avalanche 
will occur or how large it will be. 
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It has recently been shown that the energy release due to solar particle events is 
consistent with the dynamics of a SOC system [Xa06]. This was based on three analyses of 28 
years of solar proton data taken by the IMP-8 and GOES series of satellites. The first was 
rescaled range (R/S) analysis, which was used to determine that events show long-term 
correlation. The second was a demonstration of fractal properties of event sizes, which suggests 
scale invariant behavior. The third was an analysis of the number density distribution of fluence 
magnitudes, which was shown to be a power function. These are hallmark features of systems 
that exhibit self-organized criticality. 

The third of these analyses is a necessary characteristic of SOC phenomenon [Ba96], 
[Je98], [Pe02]. The number density distribution of monthly solar proton fluences for a 28-year 
period is shown in Figure 12. The ordinate represents the number of occurrences when the 
monthly fluence exceeds that shown on the abscissa. It is seen that this distribution is a straight 
line on a semi-logarithmic plot that spans about 4 orders of magnitude. The number density 
function is [Xa06] 


dN _ -29.4 
dO ~ O 


( 1 ) 


where <E> is the monthly fluence. 

In this case the density function turns out to be exactly proportional to the reciprocal of 
the fluence. Thus, the solar event data can be represented by a power function of a type 
commonly referred to as 1/f [Ba87]. It can therefore be viewed as 1/f noise, also known as 
flicker noise. It is well known that this type of noise results when the dynamics of a system is 
strongly influenced by past events. Thus, an especially compelling argument can be made that 
solar particle events are a SOC phenomenon. 
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Monthly Solar Proton Fluence (p/cm 2 -sr) 


Figure 12. Integral distribution of monthly solar proton fluences >1.15 MeV, from 1973 to 2001 

[Xa06], 


The general behavior of a SOC system is that of a non-equilibrium system driven by a 
slow continuous energy input that is released in sudden bursts with no typical size as indicated 
by the power function distribution shown in equation (1). Although research involving SOC is 
still a developing field and there is much yet to be learned about the sun’s dynamics [Lu93], 
[Bo99], [Ga03], these results strongly suggest that it is not possible to predict that a solar particle 
event of a given magnitude will occur at a given time. It is therefore a reasonable approach to 
model solar particle events as a probabilistic phenomenon. 


4.2.2 The Maximum Entropy Principle and the Distribution of Solar Proton Event 

Magnitudes 

Given that the occurrence of solar particle events is a stochastic phenomenon, it is 
important to accurately model the distribution of event magnitudes. However, in general it can 
be rather difficult to select a probability distribution for the situation where the data are limited. 
There have been a number of empirical assumptions that the event magnitudes can be 
represented by certain distributions. For example, lognormal distributions [Ki74], [Fe90] and 
power function distributions [Ga96], [Ny99] have been used. The lognormal distribution 
describes the large events well but underestimates the probability of smaller events. On the other 
hand power functions describe the smaller events well but overestimate the probability of larger 
events. This section describes a method for making arguably the best selection of a probability 
distribution for a limited set of data that is compatible with known information about the 
distribution. 
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The Maxi m um Entropy Principle was developed by E.T. Jaynes [Ja57] using the concept 
of entropy originated by Shannon [Sh49]. Jaynes showed in his studies of statistical mechanics 
that the usual statistical distributions of the theory could be derived by what became known as 
the Maximum Entropy Principle. This led Jaynes to re-interpret statistical mechanics as a form 
of statistical inference rather than a physical theory. It established the principle as a procedure 
for making an optimal selection of a probability distribution when the data are incomplete. 
Entropy is defined mathematically the same way as in statistical mechanics but for this purpose it 
is a measure of the probability distribution’s uncertainty. The principle states that the 
distribution that should be selected is the one that maximizes the entropy subject to the 
constraints imposed by available information. This choice results in the least biased distribution 
in the face of missing information. Choosing the distribution with the greatest entropy avoids the 
arbitrary introduction or assumption of information that is not available. It can therefore be 
argued that this is the best choice that can be made using the available data. 

The probability distribution’s entropy, S, is defined [Ja57], [Ka89] 

S = -\p(M)ln[p(M)]dM (2) 

where p(M) is the probability density of the random variable M. For the case of solar particle 
event fluences, M is conveniently taken as the base 1 0 logarithm of the event fluence. A series 
of mathematical constraints are imposed upon the distribution, drawing from known information. 
In this case the constraints are [Xa99]: 

a) The distribution can be normalized. 

b) The distribution has a well-defined mean. 

c) The distribution has a known lower limit in the event fluence. This may correspond to a 

detection threshold, for example. 

d) The distribution is bounded and consequently infinitely large events are not possible. 

The resulting system of equations are used along with equation (2) to find the solution 
p(M) that maximizes S. This has been worked out for many situations [Ka89] and can also be 
solved using the LaGrange multiplier technique [Tr61]. Using this procedure the following 
result for solar proton event fluences has been obtained for the solar maximum time period: 


N = N„ 


r b - 


1 -b 
^max 


(h . — (h 

t mm t ] 


-b 

max 


( 3 ) 


where N is the number of events per solar maximum year having a fluence greater than or equal 
to </>, N tot is the total number of events per solar maximum year having a fluence greater than or 
equal to <f> mm , -b is the index of the power function, and <f> ma x is the maximum event fluence. 
Equation (3) is a truncated power function in the event fluence. It behaves like a power function 
with an index of -b for (f> « <f>, max and goes smoothly to zero at the upper limit </> max - 

Figure 13 shows > 30 MeV' solar proton event data compared to the best fit to 
equation (3). The data are from the 21 solar maximum years during solar cycles 20 - 22. It is 
seen that the probability distribution derived from the maximum entropy principle describes the 
data quite well over its entire range. This strong agreement indicates that this probability 
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distribution captures the essential features of a solar proton event magnitude distribution. It is a 
power function for small event sizes and falls off rapidly for very large events. The 
interpretation of the maximum fluence parameter </> max is interesting in itself and will be discussed 
further in section 4.2.5. 



10 6 10 7 10 8 10 9 10 10 
> 30 MeV Fluence (cnr 2 ) 


Figure 13. Comparison of the maximum entropy theory result for the distribution to 3 solar 

cycles of data during solar maximum [Xa99]. 


4.2.3 Cumulative Fluence During Solar Maximum 

During a space mission the solar particle event fluence that accumulates during the solar 
maximum time period is often the dominant contribution to the total fluence. Thus, much prior 
work focuses on this period of the solar cycle. A solar cycle typically lasts about 1 1 years. A 
commonly used definition of the solar maximum period is the 7-year period that spans a starting 
point 2.5 years before and an ending point 4.5 years after a time defined by the maximum 
sunspot number in the cycle [Fe93]. The remainder of the cycle is considered solar minimum. 

Once the initial or underlying distribution of event sizes during solar maximum such as 
that shown in Figure 13 is known, it can be used to determine the accumulated fluence for a 
period of time during solar maximum. Due to the stochastic nature of the events, confidence 
level approaches are often used so that risk-cost-performance tradeoffs can be evaluated by the 
designer. The first such model was based on King’s analysis of >10 to >100 MeV protons 
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during solar cycle 20 [Ki74], [St74]. One “anomalously large” event, the well-known August 
1972 event, dominated the fluence of this cycle so the model predicts the number of such events 
expected for a given mission length at a specified confidence level. Using additional data, a 
model from JPL emerged in which Feynman et al. showed that the magnitude distribution of 
solar proton events during solar maximum is actually a continuous distribution between small 
events and the extremely large August 1972 event [Fe90], Under the assumptions that this 
underlying distribution can be approximated by a lognormal distribution and that the occurrence 
of events is a Poisson process, the JPL Model uses Monte Carlo simulations to calculate the 
cumulative fluence during a mission at a given confidence level [Fe90], [Fe93]. An example of 
this is shown in Figure 14 for > 30 MeV protons. Thus, according to this model, there is 
approximately a 10% probability of exceeding a proton fluence of 10 10 cm' 2 for a 3 -year period 
during solar maximum. This corresponds to a 90% confidence level that this fluence will not be 
exceeded. 



Fluence (crrr 2 ) 

Figure 14. JPL91 solar proton fluence model for > 30 MeV protons. The misprint of x-axis 
units has been corrected from the original reference [Fe93]. 


An underlying assumption of the JPL Model is that the number of events from year to 
year during solar maximum has such a large variation that each year during solar maximum can 
be treated the same. Other Monte Carlo based models have been developed that parameterize the 
number of events that are predicted to occur as a function of time from the beginning of the solar 
cycle. One such model, developed at ONERA, uses measured event numbers for this [Ro05]. 


1-20 




To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components 21 

and Systems (RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 

Another such model assumes the event number is directly proportional to the sunspot number 
and thus relies on knowledge of sunspot numbers during the mission time period [Ny99]. 

It has also been demonstrated that the cumulative fluence distribution during solar 
maximum is consistent with a lognormal distribution for periods of time up to at least 7 years 
[XaOO]. This was shown using the Maximum Entropy Principle, Bootstrap-like methods [Ef93] 
and by Monte Carlo simulations using the initial distribution shown in Figure 13. Thus the 
cumulative fluence distribution is known once the parameters of the lognormal distribution are 
determined. These parameters depend on the proton energy range and the mission duration. 
They have been determined from the available satellite data and well-known relations for 
Poisson processes. Figure 15 shows examples of the annual proton fluences for >1, >10 and 
>100 MeV protons plotted on lognormal probability paper. This figure is constructed so that if a 
distribution is lognormal, it will appear as a straight line. The fitted data can also be used to 
determine the lognormal parameters for different periods of time and is used in the ESP Model 
[Xa99a], 



Cumulative Probability 

Figure 15. Cumulative annual solar proton event fluences during solar maximum periods for 3 
solar cycles plotted on lognormal probability paper. The straight lines are results for the ESP 

model [XaOO]. 


Figure 16 shows a representative comparison of the models discussed above. In addition 
it shows an update of the ESP Model, called PSYCHIC [Xa04], in which the data were extended 
to cover the time period from 1966 to 2001 and the proton energy range extended to over 300 
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MeV. Results shown are for the 90% confidence level and for a mission length of two solar 
maximum years. In all cases the energy range shown corresponds to the data range on which the 
statistical models are based, i.e. no extrapolations are used. Thus, the model differences seen are 
an indicator of model uncertainties. The spectral shape for the King Model is based on the 
August 1972 event and is therefore somewhat different than the other model results. The JPL91, 
ESP, and PSYCHIC models all agree reasonably well for their common 1 to 60 MeV energy 
range. 



Figure 16. Comparison of different models of cumulative solar proton event fluence during solar 
maximum for a 2-year period and the 90% confidence level [Xa04]. 


4.2.4 Cumulative Fluence During Solar Minimum 

It has often been assumed that the solar particle event fluence during the solar minimum 
time period can be neglected. However, for missions that are planned mostly or entirely during 
solar minimum it is useful to have guidelines for solar particle event exposures, especially 
considering the current frequent use of commercial-off-the-shelf (COTS) microelectronics, 
which can exhibit rather low TID failure levels. 

Due to the relative lack of events during solar minimum, models are more difficult to 
construct for this period. However, Monte Carlo based models that parameterize the number of 
events that are predicted as a function of time throughout the solar cycle are useful for predicting 
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cumulative fluences during solar minimum. Two such models are those developed at ONERA 
[Ro05] and by Nymmik [Ny99]. In addition, recent solar minimum time periods have been 
analyzed to obtain 3 average solar proton flux levels that allow varying degrees of conservatism 
to be used [Xa04]. 


4.2.5 Extreme Value Theory and Worst Case Events 

An important consideration for spacecraft designers is the worst-case solar particle event 
that occurs during a mission. One approach is to design to a well-known large event such as that 
which occurred in October 1989 [Ty97], or a hypothetical one such as a composite of the 
February 1956 and August 1972 events [An94], Energy spectra of some of the most severe solar 
proton events during solar cycles 19-22 are shown in Figure 17. In addition, there are event 
classification schemes in which the magnitudes range from “small” to “extremely large” that can 
be helpful for design purposes [St96], [Ny96]. 



Kinetic Energy (MeV) 

Figure 17. Some of the most severe solar proton event energy spectra in solar cycles 19-22 

[Wi99], 


However, more useful information can be provided to the designer if a confidence level 
associated with the worst case event is known for a given mission length. The designer can then 
more systematically balance risk-cost-performance tradeoffs for the mission in a manner similar 
to what is done for cumulative fluences. Once the initial probability distribution such as that 
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shown in Figure 13 is determined it becomes possible to construct such a statistical model using 
extreme value theory. 

In the usual central value statistics, the distribution for a random variable is characterized 
by its mean value and a dispersion indicator such as the standard deviation. Extreme value 
statistics, pioneered by Gumbel [Gu58], focuses on the largest or smallest values taken on by the 
distribution. Thus, the “tails” of the distribution are the most significant. For the present 
applications the concern is with the largest values. An abbreviated description of a few useful 
relations from extreme value theory is given here. Further detail can be found elsewhere [Gu58], 
[An85], [Ca88]. 

Suppose that a random variable, x, is described by a probability density p( x) and 
corresponding cumulative distribution P(x). These are referred to as the “initial” distributions. If 
a number of observations, n, are made of this random variable, there will be a largest value 
within the n observations. The largest value is also a random variable and therefore has its own 
probability distribution. This is called the extreme value distribution of largest or maximum 
values. These probability distributions can be calculated exactly. The probability density is 

/max (*; n ) = n [ p ( x )X~ X P( x ) (4) 

and the cumulative distribution is 

*•«(*;») = [«*)! (5) 

An example of the characteristics of such a distribution is shown in Fig. 1 8 for ^-values 
of 10 and 100 compared to the initial distribution (n = 1), taken to be Gaussian. Note that as the 
number of observations increase the distributions become more highly peaked and skewed to the 
right. 
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Random Variable (Arbitrary Units) 

Figure 18. Extreme value distributions for n - values of 10 and 100 compared to the initial 

Gaussian distribution [Bu88], 


As n becomes large, the exact distribution of extremes may approach a limiting form 
called the asymptotic extreme value distribution. If the form of the initial distribution is not 
known but sufficient experimental data are available, the data can be used to derive the 
asymptotic extreme value distribution by graphical or other methods. For practical applications 
there are 3 asymptotic extreme value distributions of maximum values - the type I or Gumbel, 
type II and type III distributions. 

Examples of extreme value modeling of environmental phenomena such as floods, wave 
heights, earthquakes and wind speeds can be found in a number of places [Gu58], [An85], 
[Ca88]. This modeling was first applied to radiation effects problems by Vail, Burke and 
Raymond in a study of high density memories [Va83]. It has turned out to be a very useful tool 
for studying the response of large device arrays to radiation. One reason is that the array of 
devices will fail over a range of radiation exposures and it is important to determine at what 
point the first failure is likely to occur. Other radiation effects applications have been found for 
arrays of gate oxides [Va84], [Xa96], sensor arrays [Bu88], [Ma89] and EPROMs [Mc00]. 

For the application to solar particle events the interest is in the worst-case event that will 
occur over a period of T solar maximum years. Since the number of events that can occur over 
this period is variable, the expression for the extreme value distribution must take this into 
account. Assuming that event occurrence is a Poisson process [Fe93], it can be shown that the 
cumulative, worst case distribution for T solar maximum years is [Xa98a] 
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^(^;n = exp{-W„r[l-P(M)]j (6) 

where P(M) is the initial cumulative distribution, which is closely related to equation (3) [Xa99]. 

Figure 19 shows results for worst-case event fluences for mission lengths of 1, 3, 5 and 
10 solar maximum years. The ordinate represents the probability that the worst-case event 
encountered during a mission will exceed the > 30 MeV proton fluence shown on the abscissa. 
Also shown in the figure by the vertical line denoted by “Design Limit” is the maximum event 
fluence parameter, <j) max . As will be discussed next, this parameter can be used as an upper limit 
guideline. Results analogous to these have also been obtained for peak solar proton fluxes 
during events [Xa98], which are very relevant for SEE. The event fluence magnitudes are 
discussed here because of the interesting comparison that can be made with historical data to 
help validate the model. 



Figure 19. Probability model for worst-case event fluences expected during the indicated time 

periods during solar maximum [Xa99]. 


A unique feature of this model is the upper limit parameter for a solar proton event 
fluence, <j) max . For the case of > 30 MeV protons this turns out to be 1.3 x 10 10 cm' 2 . However, 
this is a fitted parameter that was determined from limited data. There must be some amount of 
uncertainty associated with the parameter. Thus, it should not be interpreted as an absolute 
upper limit. One method of estimating its uncertainty is the parametric “bootstrap” technique 
[Ef93]. This method attempts to assess the uncertainty of the parameter due to the limited nature 
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of the data. The idea is to randomly select event fluences according to the distribution given by 
equation (2) until the number of events in the distribution is simulated. The equation is then 
fitted to the simulated data, and the parameters extracted. The procedure is repeated, and each 
time the parameters have different values. After a number of simulations, the standard deviation 
of the parameter of interest can be determined. This technique showed the upper limit parameter 
plus one standard deviation equaled 3.0 x 10 10 cm' 2 [Xa99]. 

A reasonable interpretation for the upper limit fluence parameter is that it is the best 
value that can be determined for the largest possible event fluence, given limited data. It is not 
an absolute upper limit but is a practical and objectively determined guideline for use in limiting 
design costs. 

Constraints on the upper limit of solar proton event sizes can be put on models as a result 
of studies of historical-type evidence. Relatively small fluctuations of 14 C observed in tree rings 
over a long period of time [Li80] and measured radioactivity in lunar rocks brought back during 
the Apollo missions [Re97] are consistent with the upper limit parameter but are not especially 
restrictive. The strictest constraint to date comes from analysis of approximately 400 years of 
the nitrate record in polar ice cores [McOl]. The largest event reported was estimated to be 1.9 x 
10 10 cm' 2 for > 30 MeV protons. This was the Carrington event that occurred in September 
1859. Figure 20 shows a bar graph of the upper limit parameter, (j) max , for > 30 MeV protons 
including the one standard deviation uncertainty that was estimated from the parametric 
bootstrap method. This is compared with the reported value for the Carrington event. It is seen 
that these quantities are well within the uncertainties. Also shown for reference is the value for 
the October 1989 solar particle event that is commonly used as a worst-case event. 


E 

o 

Q. 

O 

O 


CD 

O 

c 

CD 

13 


C 

CD 

> 

LU 

> 

CD 


O 

CO 

A 



- — ' ‘ . ■■■■ . 

Model Maximum Oct. 1989 Carrington 


Figure 20. Comparison of the > 30 MeV solar proton event fluences of the October 1989 event, 
the 1859 Carrington event as determined from ice core analysis [McOl], and the model upper 
limit parameter plus one standard deviation shown by the error bar [Xa99]. 
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4.3 Solar Heavy Ion Models 

Solar heavy ion models are generally not as advanced as solar proton models due to the 
large number of heavy ion species, which complicates measurements of individual species. For 
microelectronics applications, solar heavy ion models are needed primarily to assess SEE. 
Astronaut exposure is also a serious concern for manned missions. 

4.3.1 Cumulative Fluences 

One quantity of interest is the average SEE rate during a mission. This means that 
models for cumulative solar heavy ion fluence must be developed to be used along with GCR 
Models. Tylka et al. used a Monte Carlo procedure similar to the JPL91 solar proton model 
[Fe93] to predict cumulative fluences for certain elements during a mission at a specified 
confidence level [Ty97a]. This was done for 2 broad energy bins each for alpha particles, for the 
CNO group, and for Fe. It is based on the University of Chicago Cosmic Ray Telescope (CRT) 
data taken between 1973 and 1996. 

The most complete model for cumulative solar heavy ion fluences is the PSYCHIC 
Model [Xa07]. Here measured alpha particle energy spectra are scaled to proton energy spectra 
based on measurements from the IMP-8 and GOES instrumentation during the time period 1973 
to 2001. The energy spectra of remaining major heavy elements - C, N, O, Ne, Mg, Si, S and Fe 
- are scaled to the alpha particle energy spectra using measurements of the Solar Isotope 
Spectrometer (SIS) onboard the ACE spacecraft over the most recent 7 year solar maximum 
period. An abundance model is used for the remaining minor heavy elements. It is based on 
measurements from the International Sun-Earth Explorer-3 (ISEE-3) spacecraft and current 
knowledge of solar photospheric abundances and processes. Results for differential fluence- 
energy spectra are shown in Figure 21 for some of the major elements and a summed spectrum 
for atomic number Z > 28. Also shown by the points in the figure are cumulative fluence results 
for alpha particles and iron for the same conditions based on the modeling effort of Tylka 
[Ty97a], 
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Energy (MeV/nucleon) 

Figure 21. Differential fluence-energy spectra for protons, alpha particles, oxygen, magnesium, 
iron and summed spectra for Z > 28 elements for a 2-year mission during solar maximum at the 
90% confidence level. Lines are spectra calculated with the PSYCHIC model [Xa07] and points 
are obtained from [Ty97a]. 

4.3.2 Worst Case Events 

In an attempt to model worst-case events, the original CREME model [Ad87] and 
subsequently the CHIME model [Ch94] scaled heavy ion abundances to protons for individual 
events. However, this assumption that individual events with the highest proton fluxes should 
also be heavy ion rich turned out to be inconsistent with subsequent data [Re99] and led to 
worst-case event models that were too conservative [Mc94]. Modifications of the original 
CREME code were made in the MACREE model [Ma95] to define a less conservative worst- 
case solar particle event. MACREE gives the option of using a model based on the measured 
proton and alpha particle spectra for the well-known October 1989 event and an abundance 
model that is 0.25 times the CREME abundances for atomic numbers, Z > 2. A model that 
originated at JPL [Cr92] characterizes the distribution of 1 to 30 MeV per nucleon alpha particle 
event fluences using a lognormal distribution in order to assign confidence levels to the event 
magnitudes. The alpha particle data are based on measurements from the IMP-8 satellite for 
solar maximum years between 1973 and 1991. For ions heavier than Z = 2 an abundance model 
is used and the fluxes are scaled to the alpha particle flux for a given confidence level [Mc94]. 
The current version of the widely used CREME code, CREME96, uses the October 1989 event 
as a worst-case scenario. It provides 3 levels of solar particle intensity [Ty97]. These are the 
“worst week”, “worst day” and “peak flux” models, which are based on proton measurements 
from the GOES-6 and -7 satellites and heavy ion measurements from the University of Chicago 
CRT on the IMP-8 satellite. The most extensive heavy ion measurements in the model are for C, 
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O and Fe ions [Ty96]. It is noteworthy that the energy spectra of these 3 elements extend out to 
roughly 1 GeV per nucleon. The remaining elemental fluxes are determined from a combination 
of measurements limited to 1 or 2 energy bins and abundance ratios. 

Comparisons to the CREME96 worst case models have been made with data taken by the 
Cosmic Radiation Environment DOsimetry (CREDO) Experiment onboard the Microelectronics 
and Photonics Test Bed (MPTB) between 2000 and 2002 [Dy02]. The data show that 3 major 
events during this time period approximately equaled the “worst day” model. An example of this 
is shown in Figure 22 for an event that occurred in November 2001. 



100 1000 10000 
Linear Energy Transfer (MeV/(g*cnr 2 )) 

Figure 22. Comparison of a solar heavy ion event that occurred in November 2001 with the 
CREME96 “worst day” model. The progression of daily intensities is indicated with the peak 
intensity occurring on day 2929 of the mission [Dy02]. 


A summary of the heavy ion space environment is shown in Figure 23 for the solar 
maximum time period. Plotted are results for integral LET spectra for 3 situations that need to 
be considered for both spacecraft design and SEE rate predictions. Results obtained from 
CREME96 are for the GCR flux and for the “worst day” solar particle event model [Ty97]. 
Results obtained from the PSYCHIC model are for cumulative solar particle event exposure at 
the 90% confidence level [Xa07]. 
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Figure 23. LET spectra obtained from CREME96 [Ty97] for the “worst day” solar particle event 
(top curve) and GCR during solar maximum (bottom curve), compared to the cumulative solar 
particle event flux at the 90% confidence level. All results assume 100 mils of aluminum 
shielding. 


5 Earth radiation belts 

5.1 Overview and background 

5.1.1 The Earth’s magnetosphere 

The Earth's magnetosphere can be seen as a natural cavity in the interplanetary medium 
in which the Earth is relatively well protected against external influences. It is compressed on the 
solar side and highly extended on the anti-solar side. In this structure, at the level of the poles, 
two horns, flaring out towards space offer the particles from the interplanetary medium a 
possibility of penetrating into the upper atmosphere. Close to Earth, the charged particles present 
in the magnetosphere can be trapped by the magnetic field and form the radiation belts. 

In the magnetosphere, the radiation belts only occupy a relatively restricted internal 
region (Figure 24). The region closest to Earth is well known and constant over time: this is the 
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upper atmosphere. The external limit, however, is poorly defined and depends on the conditions 
in the solar wind through the magnetic field as we will see later on. They therefore extend from 
the upper atmosphere (some hundreds of km) up to geostationary orbit and beyond. 



Figure 24 - The Earth's magnetosphere. 

5.1.2 The Earth magnetic field 

In the Earth's magnetosphere, the magnetic field is the sum of two terms, one of internal 
(main component) and the other of external origin. The internal magnetic field is probably due to 
the convection motion in the core of the planet; in addition to this main term, there is the 
permanent residual field of the Earth's crust. At the zero order the field can be considered to be 
dipolar. The corresponding field lines are shown in Figure 25. However the single dipolar 
approximation is not rigorous. It is then more appropriate to take an off-center and tilted dipolar 
magnetic field as approximation. This gives a dipole whose center is not at the center of the 
Earth and whose axis is not parallel to the Earth's rotation axis (Figure 26). The result of this 
geometry of the magnetic field is an anomaly, a zone in which the field is weaker. This region is 
situated at the level of Brazil, and is known as the South Atlantic Anomaly. More realistic 
models of the internal field exist, the old models such as Jensen and Cain 1962 [Je62], GSFC 
12/66, and the International Geophysical Reference Field model (IGRF) [Ca67]. These models 
consider the terms of a multipolar higher order. The potential, V, is then a development in series 
based on Gauss spherical harmonic coefficients, gH,, h!h. The general form is then given by: 

N M ( Y + V 

B=-gradV (7) with —) [gll'cosm^+hfl'smm^prcosd 

«=1 m =\ T 

It is important to note that the Earth's field is subject to long-term changes (secular 
drifts), in particular the South Atlantic Anomaly is drifting south-eastwards. At the present time, 
we note: 
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- a decrease in the intensity of 27 nT/year (0.05 % a year), 

- a drift of the axis, resulting in a westward rotation of the southern end of the dipole 
(0.014° a year) and an increase in the shift towards the West Pacific close to 3 km a year. 



Figure 25 - Dipolar magnetic field lines. 



Figure 26 - Dipolar magnetic field tilted and off-center with respect to Earth. 


An iso-contour cartography (with a constant B modulus) for a given altitude makes it 
possible to see the South Atlantic Anomaly and its drift over time (Figure 27). 
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Figure 27 - Iso-contour of the Earth's magnetic field at an altitude of 800 km in 1965 on 
the left and in 2000 on the right deduced from the IGRF model (NP for “North Pole”, SP for 
“South Pole” and SAA for “South Atlantic Anomaly”). The eastward and southward drift of the 
South Atlantic Anomaly can be seen. 

Up to a distance of some Earth radii, 5 to 6, the magnetic field is close to a dipolar field 
and the magnetosphere is more or less in revolution. Beyond that distance, the external fields 
become less and less negligible and contribute to the deformation of the dipolar internal field. 
These external fields are the sum of several components: the compression on the day side and the 
blast in the anti-solar direction, forming the tail of the magnetosphere. Several models of external 
field exist, we can mention for example the models of Tsyganenko [Ts87, Ts89] (see Figure 28), 
Tsyganenko and Stem [Ts96], Olson and Pfitzer [0177, 0188], and Alexeev et al. [A100], These 
latter - the sum of the fields transported by the solar wind and induced by the currents in the 
magnetosphere - are subject to rapid variations. The variations in the interplanetary environment 
have an impact on the magnetosphere. The variations in the speed of the solar wind (400 to 1000 
km/s) and therefore of the energy transported are, depending on the orientation of the 
interplanetary magnetic field, more or less well transferred to the magnetosphere, increasing the 
instabilities of the external magnetic field. For example the compression of the sub-solar zone 
may be sufficient to place a geostationary satellite temporarily beyond the magnetopause; 
likewise, during these geomagnetic storms, injections of high-energy particles are observed in 
the radiation belts. 

In order to understand and reproduce the dynamics of the charged particles present in the 
magnetosphere, it is common to define magnetic coordinates (Figure 29). r is the distance from 
the center of the dipole to the point under consideration, X its latitude (and 0 its colatitude: 0 = 
ti/2 - X), and cp its magnetic longitude. A field line (or force line) is defined by the Mcllwain 
parameter, L [Mc61] roughly equal (only true with a dipolar field) to the distance (expressed in 
planet radii) from the center of the planet to the intersection point of that force line with the 
magnetic equatorial plane. A point on a force line is defined by the B parameter, modulus of the 
magnetic field at the point under consideration. B and L then represent a coordinates system 
linked to the model of the magnetic field under consideration. 


1-34 




To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components 
and Systems (RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 


35 



Figure 28 - Earth's magnetic field, external component. Tsyganenko 1987 model [Ts87]. 


B 



Dipole axis 


Figure 29 - Magnetic coordinates. 

5.1.3 Charged particles motion 

All charged particles immersed in an electromagnetic field will be subject to the Lorentz 
force: f=c^aB+e) where q is the particle's charge, v its speed, B the magnetic field and E the 
electric field. If the magnetic field is very strong and the energy of the particles is great (and 
therefore their speed too) then the effect of the electric field can be ignored and the Lorentz force 
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is reduced to F«^(va.b) .Under these conditions, the movement of the high-energy particles can be 
generally broken down into three basic periodic movements. 

Gyration: a charged particle immersed in a magnetic field will have a rotation movement 
around the field line. This movement is called gyration (Figure 30). It is then possible to define 
some magnitudes relative to this movement: 

- the Larmor radius, r _ wv i where m is the relativistic mass of the particle, v± the 

z qB 

component perpendicular to the magnetic field of the particle speed, q its charge and B the 
modulus of the magnetic field; 

2 

- the relativistic magnetic moment, u = mv± = constant , 

2 B 

Bounce: If a particle has one component of its speed parallel to the magnetic field then it 
will move along the field line (Figure 31). When making any movement the particles keep their 
relativistic magnetic moment, p constant. Since the magnetic moment has to remain constant, 
the particle which moves from the equator (point where the magnetic field is weakest along the 
field line) towards the higher latitudes will see an increasingly strong magnetic field. It is 
necessary that the perpendicular component of the speed should increase in order for p to remain 
constant. This will be possible until the perpendicular speed is equal to the particle's total speed, 
the parallel speed then being null. At this particular point the particle stops, it is at its mirror 
point. A weak force due to the gradient of the magnetic field enables this particle to go 
backwards to its other mirror point situated in the other hemisphere. The particle therefore has a 
back and forth movement between its two mirror points, this is the bounce movement. 

It is possible to define the angle that the particle's speed vector must have with respect to 
the magnetic field when it crosses the equator such that its mirror point is in the upper 
atmosphere. The particle is then lost and will not be able to come back. This allows us to define a 
loss cone, that is to say if the speed vector is within the cone then the particle cannot bounce and 
will be lost. 

Drift: in order to simplify the problem, we place ourselves on the magnetic equator 
(Figure 32). Since the magnetic field of the planets has a radial gradient, the gyration cannot take 
place in a constant Larmor radius. Indeed, the magnetic field along a gyration becomes stronger 
if the particle approaches the planet, the Larmor radius is then smaller and therefore the radius of 
the trajectory's curve is also smaller. The particle will thus be able to move away from the planet, 
the magnetic field will be weaker and therefore the Larmor radius and the radius of the 
trajectory's curve will be greater. The particle therefore does not go through a simple circle but 
along a more complex trajectory. This movement breaks down into a simple gyration (circular) 
and a rotation movement around the planet: this is the drift movement. 
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Figure 3 1 - Bounce movement of a charged particle between its two mirror points 


A charged particle submitted to these three basic [No63] and periodic movements then 
moves through torus shaped surfaces around the Earth, which are commonly called drift shells 
(Figure 33). The periods associated with each of these basic movements for a 3 MeV electron at 
L=3 are respectively 2.14 10' 4 s, 0.19 s and 504 s. The disparity between the periods is very 
great, a factor of the order of 1000 should be noted between each of them going from the 
gyration movement to the drift movement. 
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Figure 32 - Drift movement for an equatorial charged particle (which does not bounce 



Figure 33 - Composition of a charged particle's three periodic movements: gyration, 
bounce and drift. The particle then follows a torus surface called a drift shell. 


5.2 Description of radiation belts 
5.2.1 Composition and topology 

The magnetic field in the vicinity of the Earth becomes such that all relativistic charged 
particles are trapped and their movement is then quasi-periodic. These special conditions are thus 
favorable to the accumulation of high-energy charged particles in certain regions of space which 
creates the radiation belts. Given the trajectories of the particles the radiation belts have a 
toroidal shape which surrounds the Earth. The Earth's atmosphere is the lower limit of the 
radiation belts since it causes the loss of all the trapped particles. The upper limit, however, is 
less clear and is defined by the minimum intensity in the presence of disturbances of the 
magnetic field such that the particles are always trapped. 
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Discovered during the first space missions by J. Van Allen, the particles trapped in the 
radiation belts (or Van Allen belts) are essentially protons and electrons. The energy ranges 
commonly encountered go from some keV up to some tens or even hundreds of MeV. Table 3 
summarizes the properties of the Earth's radiation belts. 



Particle 

Energy 

Extension (Earth 
radii) 

L Q W-h 

e- 

1 keV-7 MeV 

1-10 


p+ 

1 keV-300 MeV 

1-7 


Table 3 - Characteristics of the Earth's radiation belts. 

A view of the radiation belts is given in the following figures. A single maximum is 
observed for the proton belt (Figure 34) for a value of L that depends on the energy (L=1.7 for 10 
MeV protons); the flux is very stable there and the maximum energies can reach between some 
MeV and some hundreds of MeV depending on the position. 

The electron belt is more complex (Figure 35) and has two maximums respectively 
corresponding to the internal and external zones: 

- the first one centered on L = 1.4 extends up to L = 2.8; the electron populations are 
relatively stable there and can reach maximum energy levels of the order of 10 or even 30 MeV; 

- the second one, centered on L = 5, extends from L = 2.8toL=10; the electron flows 
there are much more variable and the energy levels can be as high as 7 MeV. 
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Figure 34- Proton radiation belt 
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Figure 35- Electron radiation belt. 


At zero order the radiation belts can be considered to be symmetrical in longitude in a 
region going from the Earth's surface up to the geostationary orbit, i.e. as long as the magnetic 
field is not too different from a dipole. However, as demonstrated in section 5.1, at high 
altitudes, the field differs from a dipole and the belts are no longer axisymmetric. At the level of 
the geostationary orbit, the fluxes of high-energy particles (electrons between 100 keV and some 
MeV and protons between 100 keV and 1 MeV) then have a maximum on the day side and a 
minimum on the night side. We then speak of day-night asymmetry (Figure 36) due to the 
topology of the magnetic field (external). 
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Figure 36- Electron fluxes with an energy level between 315 and 500 keV in 
geostationary orbit measured by Los Alamos National laboratory (NM- USA) satellite 1989-046. 
The dotted lines represent the passage of the satellite at local midday. 

5.2.2 The South Atlantic Anomaly 

Since the Earth's dipole is tilted and off-centered by 500 km towards the West Pacific, the 
radiation belt (protons and electrons) goes down to a low altitude over the South Atlantic, the 
populations of charged particles being attached to the magnetic field. A satellite in low orbit 
(LEO : Low Earth Orbit) will thus only be exposed to radiation on certain fractions of the orbit 
as far as the trapped particles are concerned when passing through (Figure 37) : 

- the polar horns (electrons below 1000 km, electrons and protons above that altitude), 

- the South Atlantic Anomaly (protons and electrons at all altitudes). 

As can be seen in Figure 38 the position of Kourou, close to the SAA, means that the 
launcher's trajectory passes through a zone with a great flux of energetic trapped protons when 
being injected into geostationary transfer orbit. This must be taken into account when designing 
the on-board electronics which may be sensitive to the singular effects induced by protons. 


1-41 



To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components 
and Systems (RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 


42 



Figure 37 - Environment in low orbit 
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Figure 38 - Iso-flux curves for 9.4 MeV protons (top) and for 460 keV electrons (bottom) 
measured by the ICARE detector on the Argentinean SAC-C satellite at an altitude of 710 km. 
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5.3 Dynamics of the radiation belts 


Given the measurements of the trapped particles in the Earth's environment, it is now 
certain that a static view of the radiation belts is obsolete. The American CRRES satellite in the 
early 1990s clearly evidenced the extreme dynamics of the trapped electrons and protons: As 
stated earlier, the radiation belts are linked to the existence of the Earth's magnetic field and the 
populations of particles are the result of an equilibrium between: 

- the sources, injections from the tail of the magnetosphere and creations by nuclear 
reactions between atoms in the upper atmosphere and energetic ions (solar or cosmic), 

- the losses by precipitation in the upper atmosphere or by charge exchange with the 
atoms and molecules from the exosphere (extended atmosphere). 

These various terms can vary over time and are highly dependent on the magnetic field; 
any transient disturbances and time drifts of the magnetic field result in rapid fluctuations 
(magnetic storms) and long-term variations of the fluxes in the belt. 

It is therefore judicious to look at the dynamics of belts at different time scales as a 
function of the effect (of the degradation) to be studied. If you are only interested in the 
cumulative effects such as the dose effects then variations on the scale of the minute (sub-storm) 
or of the week (storm) serve no purpose. However the variations on the solar cycle scale will be 
fundamental. On the contrary, in the case of studies on charging environments, the time scales go 
from some hours to some days, and the dynamics of the belts must be described on the scale of 
the sub-storm (surface charge) or of the magnetic storm (surface and internal charge). 

5.3.1 Dynamics on the scale of the solar cycle 

Protons 

The proton radiation belt (high-energy component >10 MeV) varies slowly as a function 
of the solar cycle [Hu98]. The flux levels are roughly at their highest when the solar cycle is at 
its lowest and vice versa. This is the result of two physical processes that condition the dynamics 
of the protons, the absorption of the protons by the upper atmosphere on the one hand and the 
modulation of the CRAND source (Cosmic Ray Albedo Neutron Decay) on the other hand. 
When the solar cycle is at its maximum the upper atmosphere is heated up and the densities at 
constant altitude increase. It can then be understood that the losses of trapped protons induced by 
the charge exchange increase. However, when the solar cycle is at its maximum, the fluxes of 
cosmic radiation fall due to the intense solar activity and the source is reduced. The balance is 
shown in Figure 39. 

Another important characteristic of low altitude proton fluxes is an overall decrease in the 
fluxes from one cycle to the next. The comparison of the fluxes from one solar minimum to the 
next shown in Figure 39 very clearly illustrates this slow variation. It is due to the secular drift of 
the Earth's magnetic field. At the present time this poses problems for forecasting the low 
altitude proton fluxes from one cycle to the next. 
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Figure 39 - Changes in the proton fluxes at low altitudes (bottom), in the cosmic radiation 
(middle) and atmospheric densities (top) as a function of the solar cycle. 

Electrons 

The variations in the electron belt are above all known in the external zone, the 
geostationary orbit being particularly well documented [Ba86, Re98]. In geostationary orbit 
(Figure 40), the fluxes of electrons are at their lowest when the solar cycle is at its highest and 
are at their highest three or four years after the top of the cycle (just before the solar cycle is at its 
lowest). This modulation as a function of the radio-solar flow (FI 0.7) at 10.7 cm shows that the 
amplitude is all the greater if we examine the high energies (MeV and above). However at low 
energy levels (some hundreds of keV) the modulation is virtually inexistent. The strong fluxes of 
electrons are linked to the presence of coronal holes on the surface of the sun which in turn 
involve intense and long-duration magnetic storms at the level of the Earth's magnetosphere (see 
next paragraph). 
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. Figure 40 - Electron fluxes at geostationary orbit as a function of the solar cycle. 


Even if the charge phenomena on satellites are induced by instantaneous fluxes of 
electrons, it is nevertheless possible to use this curve to define the unfavorable periods that lead 
to surface or internal charging. Surface charging can appear at any moment during the solar 
cycle, since the low-energy electrons involved in this process are not modulated by the solar 
activity. As for the internal charge, it will preferably appear some years after the maximum of the 
solar cycle when the coronal holes have an influence on the Earth's environment. 


5.3.2 Dynamics on the scale of the magnetic storm 

Proton 

The low-energy protons (some tens to some hundreds of keV) are very sensitive to 
magnetic storms. The fluxes of particles therefore follow the Earth's magnetic activity in a region 
going from L=2 to L=6 with time scales going from a minute to several hours. A view from the 
CRRES satellite (MEB instrument) makes it possible to view 14 months' dynamics of the belt of 
62 keV protons (Figure 41) in various different regions (the satellite crosses the magnetic field 
lines close to the equator) [Fr95, Bo98]. 

At higher energy levels (several tens of MeV) the belt is generally very stable but major 
events can dramatically change the flux levels in intermediate regions of the radiation belts 
[Gu96]. In fact, if there is a solar flare in progress when a very intense magnetic storm is 
building up then the solar particles can be trapped and thus significantly increase the fluxes of 
trapped particles in a region between L=2 and 2.5 [Va99]. The example of the event in March 
1991 is striking (Figure 42). 
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Figure 41 - Flux of 62 keV protons measured by the MEB detector on CRRES for 14 
months. * 
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Figure 42 - Flux of 36.3 MeV protons measured by the PROTEL detector on CRRES for 
14 months (march 1991 event shown by red arrow). 
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Electrons 

To better understand the environmental conditions that lead to charging phenomena on 
satellites the dynamics of the electrons will have to be controlled during major magnetic storms. 
According to the observations made in geostationary orbit the low-energy electrons (which 
induce surface charges) appear right from the first instants of the disturbance whereas the higher- 
energy electrons (which induce the internal charge) are detected some days after the beginning of 
the event (Figure 43). 



1994-084 
Elec. 50.-75 keV 


1994-084 

Elec. 1100.-1500. keV 


omni 
Kp Index 


00:00 18:00 12:00 06:00 00:00 Time (UT) 
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Figure 43 - Flux of 50-75 keV (top), 1.1-1. 5 MeV (middle) electrons in geostationary 
orbit and magnetic activity Kp (bottom) as a function of time (Peak intensity of the storm shown 
by red arrow). 


A distinction can also be made between various different classes of events [McOlb] 
(Figure 44) : (1) storms where the magnetic activity index Kp is high (up to 8 or 9) but which do 
not last long (less than 1 day) and which do not produce any or only a few high-energy electrons 
in geostationary orbit (Figure 44 on the left) and (2) the storms where the magnetic activity index 
Kp is moderate (up to 6 or 7) but which last longer (several days) and which produce large 
quantities of high-energy electrons in geostationary orbit (Figure 44 in the middle). The most 
surprising thing is that in the second case the storms where the magnetic activity index oscillates 
between 2 and 4 for several days produce nearly as many high-energy electrons in geostationary 
orbit as a more violent storm with a Kp of 6 (Figure 44 on the right). 
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Figure 44 - Comparison of the consequences of three magnetic storms on the fluxes of 
1.1 -1.5 MeV electrons in geostationary orbit for three different levels of activity 

In order to understand the electron acceleration phenomena we must examine in detail the 
physical processes that affect the dynamics of these particles. The low-energy electrons are first 
of all transported from the tail of the magnetosphere towards the inside of the radiation belts by 
an increase in the radial diffusion at the beginning of the storm. As they come closer to Earth, 
they see an increasingly strong magnetic field and, under the effect of the Lorentz force, they 
drift around the Earth firstly in the night-morning sector. In the vicinity of the plasmapause 
(~L=4), these electrons are not only submitted to radial diffusion but also to the wave-particle 
interaction. The combination of the two leads to a slow but continuous acceleration of the 
electrons. A large proportion of these electrons (now at high energy levels) will diffuse radially 
towards the internal zone of the radiation belts. There, they will be globally lost by precipitation 
in the loss cone due to the wave - particle interaction which is then preponderant. However a 
small proportion of the high-energy electrons will be able to diffuse from the plasmapause 
towards the exterior of the belt and finally reach the geostationary orbit. This acceleration 
mechanism makes it possible to explain, in particular, the delay between the low-energy 
electrons and the high-energy electrons in geostationary orbit. 

Events such as these are very frequent, the CRRES period makes it possible to view the 
dynamics of the electrons for 14 months still with the major storm of March 1991 (Figure 45). 
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Figure 45 - Flux of 1.6 MeV electrons measured by the HEEF detector on CRRES for 14 

months 


5.3.3 Extreme events in the Earth electron belts 

It is clear that there are many magnetic storm affecting the Earth radiation belts and 
some of them can be very extreme. Of course they are likely to happen rarely but can affect any 
space system dramatically. One full solar cycle is represented in Figure 46 to appreciate the 
occurrence of such strong storms where extreme electron events recorded along two different 
orbits are highlighted. The sunspot number is plotted on the top panel to identify the phase of the 
solar cycle. The middle panel shows 5.35 MeV electron fluxes measured at LEO onboard 
NPOES-15 spacecraft (800 Km- 98°) and the bottom one shows 5.5-7. 1 MeV electron fluxes 
measured along Polar orbit (HEO). Over this long time period extreme events are only recorded 
during the declining phase of the solar cycle at LEO, this was during the large July- August 2004 
and November 2004 storms. During that time very energetic electron were produced. Usually 
environment at LEO is assumed to be mainly a proton environment and any anomaly, like single 
event transient (SET) is assumed to be induced by definition by a single particle. In such a case 
what is called a SET could just be a transient anomaly induced by an internal discharge. On the 
other hand along Polar orbit the situation is different. Along this time period the extreme events 
(regarding 5.5 - 7.1 MeV electrons) are seen in August-September 1998 during the rise of the 
solar cycle. Also the location of maximum flux values are at larger L shells compared to the 
August 2004 events. As a conclusion, conditions under which an extreme electron environment 
built up at LEO can be very different to the ones at HEO orbits. In other words, when an 
"extreme event" is seen along a given orbit does not mean it is extreme for all orbits. Non 
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linearity in particle dynamics makes it difficult to extrapolate any single, local measurements to 
all locations in radiation belts. 
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Figure 46 -Top panel: Sunspot number, middle panel: 5.35 MeV electrons measured at 
LEO onboard NPOES-15 and bottom panel: 5.6-7. 1 MeV electrons measured at HEO onboard 
Polar. 


5.3.4 Extreme events in the Earth proton belts 

For proton radiation belts, extreme event results from a combination of a solar flare 
(which is from the radiation belt point of view a source term) and a large magnetic storm (in this 
case usually due to a coronal mass ejection). An example is given in Figure 47. The top panel 
shows 9-15 MeV solar protons measured by GOES-08 spacecraft at geosynchronous, middle 
panel shows in a L versus time map 9.65-1 1.35 MeV protons measured at LEO onboard SAC-C 
and the bottom shows the magnetic activity index Kp. On the middle panel proton flares are 
clearly observed at L values greater than 4 (Note that L of 4 highlighted by the horizontal dashed 
line show the average magnetospheric shielding) whereas trapped particles in the proton 
radiation belts are encountered at L below 3. On march the 31 st , a large coronal mass ejection 
impacted the Earth magnetosphere and induced a large magnetic storm. The net effect on protons 
was to lower down the magnetospheric shielding leading access to low L shells to solar protons. 
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Then a large amount of these solar particles have been trapped while the magnetospheric 
shielding was low. The trapped proton enhancement is clearly seen in the L=2 range. 
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Figure 47 -Top panel: Solar protons observed by GOES-8 at GEO, middle panel: 9.65- 
1 1.35 MeV protons measured at LEO onboard SAC-C and bottom panel: Kp index. 

Once new protons are trapped in the belts they can remain there for months. Fluxes 
can decrease slowly because particles are lost by friction processes and/or charge exchange or 
can be lost suddenly because a new large magnetic storm occurs with no source term from solar 
flares (see Figure 48). 
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Figure 48 -9.65-1 1.35 MeV protons measured at LEO onboard SAC-C. 

In other words, some extreme events combined with a solar flare can lead to large 
flux enhancements whereas others not being combined with a flare can lead to large flux 
decreases. It makes such events difficult to predict because the net results before and after the 
extreme event depends on how the storm and the flares are synchronize at the Earth. 

5.4 Static models 


5.4.1 NASA's AP8 and AE8 models 

Numerous measurements performed between 1958 and 1978 have made it possible to 
familiarize ourselves with the fluxes and energies of the electrons and protons trapped in the 
Earth's radiation belts. Empirical models have been derived from these measurements, giving a 
general but static view of the belts (Figure 49). The most recent ones, developed by the NASA in 
the 1970s, AP8 ("Aerospace Corporation Proton version 8") and AE8 ("Aerospace Corporation 
Electron version 8"), give proton and electron spectrums at the solar minimum (AE8 MIN and 
AP8 MIN) and maximum (AE8 MAX and AP8 MAX) at all geomagnetic coordinate points (B, 
L) in the magnetosphere. Their validity range extends respectively from L= 1.15 to 6.6 for 
protons with an energy level comprised between 100 keV and 400 MeV and from L=1.2 to 11 
for electrons with an energy level comprised between 40 keV and 7 MeV [Sa76, Ve91, Fu97]. 
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Figure 49 - Omnidirectional integrated fluxes of protons with an energy level higher than 
10 MeV (left) and of electrons with an energy level higher than 1 MeV (right). 


These models, however, are now obsolete even if they remain a reference for all 
industrial companies working in the space sector [Da96]. 

First, the magnetic field has drifted and the South Atlantic Anomaly is now further to the 
East and South than it was in the 1970s. This problem, however, is not necessarily a limitation. 
When the fluxes of particles along a satellite orbit are assessed with a view to calculating the 
total dose for a mission, an average must be calculated for a large number of orbits. The exact 
position of the Anomaly is then no longer of interest. However, if for a certain mission, an 
Ariane launch from Kourou for example, the exact position of the South Atlantic Anomaly is 
required, it would be judicious to use the model of the magnetic field from that time and then 
make the latitude and longitude transformations to take into account the Anomaly's drift. It has 
also been demonstrated that the fluxes of protons are underestimated (by a factor of 1.6 to 2) by 
the AP8 models for altitudes comprised between 300 and 500 km, i.e. in the vicinity of the cutoff 
induced by the atmosphere. This is partly due to the interpolation technique. Finally, the East- 
West asymmetry of the particles measured at those altitudes is not reproduced by the AP8 
models. 

The AE8 model also has its shortfalls. At low altitude, at the level of the internal belt, the 
maximum energy given by the model is 5 MeV whereas the CRRES satellite has measured 
electrons with an energy level of 30 MeV. In the outer belt, however, the flows are overestimated 
by a factor of at least 3 at high energies. Likewise in geostationary orbit, the fluxes predicted by 
the AE8 MIN and MAX models are identical, the variations due to the solar cycle are not 
correctly reproduced here. 


5.4.2 The NOAAPRO model 

The NOAAPRO (NO A A proton) model developed by S. Huston for the NASA is the first 
model that takes into account the variations of the fluxes of protons trapped in the radiation belts 
at low altitude as a function of the solar activity. It is based on the MEPED measurements 
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performed by the NOAA-TIROS satellites. It determines the mean omnidirectional integrated 
fluxes of protons with energy levels higher than 16, 30 and 80 MeV at an altitude of 800 km as a 
function of the date and of the radio-solar flow at 10.7 cm, F10.7 [Hu98]. This model is of course 
still very limited from the energy range and altitude viewpoint, however for a polar orbit at 800 
km it represents a great step forward with respect to AP8 and is therefore positioned as a good 
base for a future low-altitude proton model. 


5.4.3 The IGE2006 model 

The IGE2006 model (“International Geosynchronous Electron” previously called POLE 
“Particle Onera Lanl Environment”) developed by ONERA-DESP in cooperation with the Los 
Alamos National Laboratory (LANL) and JAXA is the first model that takes into account 
variations of the fluxes of electrons trapped in the radiation belts in geostationary orbit of a 
function of the solar activity [Bo03]. It is based on the measurements made by the LANL's and 
JAXA’s geostationary satellites. It determines the mean omnidirectional differential fluxes of 
electrons with an energy level comprised between 1 keV and 5.5 MeV as a function of the solar 
activity. Once again, this model is limited to a single orbit but, here too, it represents a good step 
forward with respect to AE8 and a first step to develop a mean model of the electron 
environment in the external belt. 

5.5 Dynamic models 

5.5.1 The CRRESPRO and CRRESELE models 

Subsequent to the CRRES mission in 1990 and 1991 (14 months in all) empirical and 
dynamic models of the radiation belts saw the light of day. The CRRESPRO (CRRES proton) 
model provides fluences of trapped protons with an energy level comprised between 1 and 1 00 
MeV for values of L between 1.15 and 5.5 (Meffert and Gussenhoven, 1994). It is based on the 
measurements provided by the PROTEL telescope. It calculates the fluences of protons averaged 
over an orbit for two states, one so-called active and another quiet state (CRRESPRO ACTIVE 
and CRRESPRO QUIET). This model reflects, in particular, the creation of a second belt of 
protons subsequent to a solar flare synchronized with a major geomagnetic storm. As for the 
CRRESELE (CRRES electron) model it provides the fluxes of electrons with an energy level 
comprised between 700 keV and 5 MeV for L values of 2.5 to 6.5 and six levels of magnetic 
activity, plus a worst case and a mean state [Br95, Br92]. It is based on the measurements given 
by the HEEF detector. The inputs are only Ap (geomagnetic activity index) averaged over 15 
days. This model makes it possible to take into account the creation of a third belt and the effects 
of recurrent storms at 27 days. The major problem of these models is their representativeness 
given the short period of measurement acquisition (14 months during a solar maximum). 


5.5.2 The ESA SEE1 model 

This model is also deduced from the measurements provided by the CRRES satellite, and 
more precisely by the MEA detector. This model provides fluxes of electrons with energy levels 
higher than 100 keV as a function of the magnetic activity index Kp [Va96]. It is associated with 
a neuronal network which predicts the fluxes on the basis of the magnetic activity index Kp 
throughout the solar cycle. This model still poses the problem of the representativeness of the 
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CRRES measurements (14 months as opposed to the 1 1 years of the solar cycle). Furthermore, 
the fluxes of electrons with energy levels higher than 1.5 MeV are extrapolated and are 
consequently overestimated. 


5.5.3 The Salammbo models 

Since the 1990s, ONERA/DESP has been developing physical models of the proton and 
electron radiation belts, called the Salammbo codes. At the present time, these codes represent a 
family of models (Salammbo 4D, 3D and 2D) which provide a more or less well-refined 
description of the belts as a function of what one wants to reproduce and of the desired resolution 
of the result [Bo96, Bo97, Bo98, Va99]. These models make it possible to describe the dynamics 
of the proton and electron belts with energy levels of 10 keV-300 MeV and 10 keV-10 MeV 
respectively in the region going from L=1 to 7 with a time resolution of between one minute and 
several hours. The inputs to these models are the magnetic activity indices Kp and Dst and a 
boundary condition deduced from geostationary measurements. At present these models make it 
possible to understand the dynamics of the charged particles trapped subsequent to magnetic 
storms of variable intensities. Notably, the creation of the second proton belt seen by CRRES has 
been reproduced, as have the effects of long and short magnetic storms on the external belts of 
electrons. In the future, this model should make it possible to define the conditions required to 
obtain a worst case for radiation belts. It also offers the possibility of validating, or not, the 
representativeness of the measurements and even of extrapolating measurements over time. 
Finally another application of the model is to interpolate and extrapolate the measurements in 
order to reconstruct a complete and dynamic cartography of the radiation belts. 

5.6 Discussion relative to the various models 

All the models that have been developed to date are based on various different data bases 
compiled at different times. It is clear that the representativeness of the measurements is a major 
problem as far as the development of models is concerned. It can be noted, in particular, that 
AP8 and AE8 are based on non-continuous measurements acquired during a small solar cycle, 
the NOAAPRO and POLE models result from data collected during solar cycles of the same 
level, and the CRRESPRO, CRRESELE and ESA SEE1 models are based on a data base that is 
extremely limited over time and only during a solar maximum (Figure 50). 
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Figure 50 - Coverage of the measurements used to create the various models of radiation 

belts. 

Another problem is the L, B and energy coverage of each of these models. The most 
complete at the present time are still AP8 and AE8. All the others only partially cover this three- 
dimensional space. A comparison is given in Figure 51. 



Energy (MeV) Energy (MeV) 

Figure 51 - Validity domain of the radiation belt models 

From the engineering model point of view, the AP8 and AE8 models are the ones that 
cover the largest domain. Updating works are in progress, they have given rise to the creation of 
new models which are still limited but which offer new bases for developing a complete model 
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of the radiation belts. Another shortfall at the present time is the definition of the worst case 
environment for assessing the transient effects that can appear on satellites. 


6 Sensitivity of orbits to the radiation 


It is clear that, given the distribution of the high-energy charged particles in the radiation 
belts and the magnetospheric shielding to protect against solar or cosmic particles, the 
environment close to satellites is highly dependent on the orbit. Here we propose to give a rapid 
overview of this environment based on the AP8, AE8 and JPL91 models for low orbits, 800 km 
98° and 30°, 1400 km circular, geostationary orbit and 20000 km 55°. 

A projection of each of these orbits in a meridian plane is shown in Figure 52 in order to 
appreciate the belt regions passed through by each of the satellites. It can immediately be seen 
that the high-altitude orbits are not subjected to the fluxes of high-energy trapped protons. 
However, these orbits are constraining from the relativistic electron viewpoint. 



Figure 52 - Projection of the various orbits in a meridian plane with, in the background, 
omnidirectional integrated fluxes of protons (AP8 MIN, energie E>10 MeV) on the left and of 
electrons (AE8 MIN, E>500 keV) on the right. The red arrows indicate the solar or cosmic 
protons with the field line where they are stopped. 


If we look in greater detail at the differences from one orbit to another (Figure 53) it can 
be seen, from the point of view of the protons, that the high orbits (35000 and 20000 km) are not 
very sensitive to the trapped protons but, on the other hand, they are not protected by the 
magnetospheric shielding, hence the presence of solar protons. In low orbit the situation is quite 
different, the trapped protons are not negligible whatever their level of energy. However, if a low 
orbit is greatly tilted then the presence of solar protons will be great, but there are virtually none 
for slightly tilted orbits, which are naturally protected by the magnetospheric shielding. As for 
the electrons, as it has been stated earlier, the high orbits are greatly exposed. 
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Figure 53 - Omnidirectional differential flux spectrums for trapped protons (AP8 min) on 
the left, for trapped electrons (AE8 max) in the center and for solar protons (Feynman 80 %) on 
the right. 


7 Future Challenges 

A number of future challenges exist for space radiation environment models if they are to 
continue to help produce reliable, cost-effective spacecraft designs and have utility in 
implementing new space technologies. There should be a goal to produce more dynamical and 
more physical models of the environment. Such results should allow more accurate projections 
for future missions. Galactic cosmic ray models are closely tied to solar activity levels, which 
modulate the fluxes of these energetic ions. Challenges for these models are to incorporate an 
improved description of the solar modulation potential and to develop cosmic ray transport 
models that account for relevant astrophysical processes. Solar particle events demonstrate a 
strongly statistical character. A major challenge for these models is to develop a description of 
the energy storage and release processes in the solar structure. This could lead to a more detailed 
probabilistic model of the cyclical dependence of event frequencies and magnitudes. 

Developing and implementing a strategy to deal with the radiation environment for 
manned and robotic space missions is critical for new interplanetary exploration initiatives. 
Getting astronauts safely to Mars and back will involve unprecedented strategies. For example, 
the lack of predictability of solar particle events indicates a potential strategy of establishing a 
measurement system in the inner heliosphere for the early detection and warning of events. 
Once an event is detected, accurate predictions must be made of the transport process to Earth, 
Mars and possibly beyond so that properties such as time of arrival, duration, intensity and 
energy spectra can be transmitted well ahead of the arrival time. The current GCR models 
depend on knowing solar activity levels in order to predict GCR fluxes. Thus, the lack of an 
established method for predicting future solar cycle activity is a concern in the planning for new 
exploration initiatives. Higher than expected GCR fluxes are a serious problem for long-term 
manned missions because they are difficult to shield against. 

For trapped particle radiations, efforts are being made to develop particle maps for various 
climatological conditions that occur throughout the solar cycle for the full range of particle 
energies and geomagnetic coordinates covered by the AP-8 and AE-8 models. Ultimately, the 
goal is to develop an accurate description of the source and loss mechanisms of trapped particles, 
including the influence that magnetic storms have on the particle populations. 


1-58 





To be presented by Mike Xapsos at 9th European Conference Radiation and Its Effects on Components 59 

and Systems (RADECS07) - Short Course Session, Monday, September 10, 2007 - Deauville, France. 

Thus, although there has been recent progress in modeling the space radiation environment, 
there are many future challenges that remain. 
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